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A КАО 


Е са formilaiion or the bed shear stress coef- 
вене inside The surf zone is derived using The concept 
Barton stress. A truncated Баутетап р.а. ғ. is used 
То describe the wave field inside the surf zone and provides 
The input to calculate the variation of wave energy and long- 
Shore current as a function of wave hetght, water depth and 
distance to shore. Тһе wave set-up is approximated using a 
sinusoidal wave solution. Field measurements of longshore 
current and waves within the surf zone are used to calculate 
the bed shear stress coefficient. The data consist of 647 
Мени ро (пто selecied from LEO program and 62 data points 
meom ingle (1966) observations, all taken along the Southern 
ПИК ог та oss i ШШ es reauency distributions and statistics 
are calculated for the bed shear stress coefficient. A mean 
bed shear stress coefficient to two significant decimal 


ШІсСсес 15 found то Бе 0.01. 
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REN TRODUCT ION 


STATEMENT 0F THE PROBLEM 

It is known that when sea waves or swell approach a 
Ки с о scoasıline at am Oblique angle a mean current is 
generated parallel to the shoreline, see Figure 1. Such 
longshore currents are of prime importance for bofh 
eoastal engineering and for aiding in the strategic planning 
Of Naval inshore warfare operaTions. 

An accepted theory of longshore currents on piane 
beaches is developed in terms of the momentum flux due to 
The waves directed down coast being balanced by the shear 
stress associated with the mean flow. The formulation of 
the bed shear stress requires the specification of a bed 
shear stress coefficient. The purpose of this thesis is 
the determination of the bed shear stress coefficient to be 
Used in The longshore current formulas. 

The study will also help in the analysis of sediment 
transport. The shear stress does work on the bottom in 
moving sediments. Several authors have formulated sediment 
transport in terms of the bed shear stress which in turn 


Beqguires an appropriate bed shear stress coefficient. 


БН HISTORICAL REVIEW 
ааа Осипа 1952) using the momentum approach for 


ШО баг сту сп от longshore current by Putnam, Munk апа 


Ш 
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EE nc i. Definition of Longshore Current Variables 
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ЕКОО О псе nar in order то fit theory with obser- 
vations, the bed shear stress coefficient must be permitted 
пи пе longshore velocity over a wide range of 

3 1/2 orders of magnitude. 

Өш Шс О 95408 Sind that the spectral limitations 
of wave growth, under the action of steady wind in shallow 
water with a typical sandy bottom, suggested a value for The 
Mate on coefficient of between 0.0l and 0.02. Also, he 
found that the observed damping of swell propagating over a 
smooth, leveled, impermeable sea bed was consistent with a 
Ка Ше от the coefficient of between 0.054 and 0.097. 

ЕЕ ЕЕ sl Using the concept of radiation 
ess, developed a relationship for prediction of the 
theoretical maximum longshore current just inside the break- 
and Broposed a friction coefficient of the order of 0.01. 
memconcitided, on rhe basis of The finding of Bretschneider 
0051). Ргапаті (1952) and Nikuradse's experiment with 
roughened pipes, that is was not "...unreasonable to expect 
ate lon cOetficient of the order of 0.01." 

БЕРІ |! (ас акеп (гот Sonu (1975); it summarizes some 
Maltes of the friction coefficient proposed by various inves- 
Tigators. The values reported were obtained from measure- 
ШОО ЫШ outside The surf zone or from laboratory experiments. 
It can be seen from this table that the range of values is 


МНО (гле (| у wide and The Test conditions varied. 
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IE IE е rieron Coefficients 


Proposed by Various 


Wave Wave 
ШЕГСТ {ton derart Period 
Coefficient (meters) (sec) 
ШЕСІ ОЕР AO | thar y 


0.030-0.089 


0.020-0.040 


0.01 


0.01 


9505 


0.05 


0.09 


-0.40 


. 


со | 2 


2070 


0.25-0051 —2,8825.96 


2 10 

5.4* 8.4% 
В 0.88-2.58 
OO 


Maz dr 9.1-15.2 


ISO 724-1252 


u Zum) оо 
U: nearbottom velocity 
n: ripple height 


Note: *Equivalent values at 10 m depth 


duo 


Test 
Cond lt ons 


Shallow water 
steady state 
wave generation 


Gul i of Mexico; 
depths 3.4-5.2 m 
slope 0.00035- 
0.00-41 


Сала авап: 
depths 2.25- 
от 

slope 0.018 


Oscillating water 
channel, turbu- 
lent boundary 
layer 


Wave flume, 
laminar boundary 
layer 


ео iB ame 
depths 13-10 m 
slope 0.0060 


Takahama, Japan; 
depths 10-7 т 
slope 0,0057 


Wave flume study; 
derived from 
energy dissipa- 
TION In Sand 
ripple vortices 


Investigators 


Authors 


Bretschneider 
(1954,а) 


Bretschneider 
(1954,b) 


Kishi (1954) 


Jonsson (1966) 


lwagaki and 
Tsuchi ya( 1966) 


Iwagaki and 
Kakinuma( 1966) 


Імадакі апа 
Kakinuma( 1966) 


Tunstall and 
Inman (1975) 





POR ECT VES Or THE STUDY 

The bed shear stress coefficients previously determined 
are based on a very limited set of field data or on labora- 
Тогу studies which used as a model simple sinusoidal waves 
which are not typical of the randomness found in nature. 
mre cop jective Of This study is fo analyze existing sets of 
field observations obtained in the surf zone and by using 
The best available theory attempt to determine a reasonable 
value of the bed shear stress coefficient. For this purpose 
a fairly large data set obtained for the Channel Island 
таста! Environment Observations (LEQ) Program was used as 
“ас а вет от оосегуатіопе by ingie (1966) taken at various 
ea Tons along the Southern California. coast. It is 
expected that the data obtained and the theory applied will 
ultimately contribute to the establishment of a reasonable 
value of the bed shear stress coefficient and to a more 
шесцгате prediction of the longshore current velocity across 


mime surf zone. 





BE AE ORY 


m INTRODUCTION 

Several models have been proposed for the distribution 
of the longshore current velocity across the surf zone on a 
plane sloping beach. The solution,which uses pure sinusoids 
to describe the waves and no lateral shear stress, gives a 
БЕУ distribution which IS triangular shaped with both 
a peak velocity and a discontinuity at the breaker point 
E (Bowen, 1969), Thornton (1969) and Longuet-Higgins (1970) |. 
This is unreasonable since there are no discontinuities in 
nature. A second model including lateral shear stress tends 
TO smooth out the discontinuity at the breaker line and pro- 
duces a smoother velocity distribution with the maximum 
ИСЕ ТУ occurring closer то shore. However, no criterion 
to predict an optimum lateral shear stress coefficient is 
as yet available. This introduces an added complexity To 
the problem. 

A random-sea model developed by Collins (1972) circum- 
Eus the difficulty of the lateral shear stress coefficients 
Ensallows the statistical input of the sea state as de- 
Eu bed by a Rayleigh distribution. Figure 2 compares the 
ВЕСТ ту distribution resulting from the various models: the 
non lateral stress model, the lateral shear stress mode! for 


EMcoestiicient equal to 0.4 (Longuet-Higgins, 1970), and the 
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random-sea model. The velocities are referenced to the peak 
ПРИ а поље поп lateral shear stress model, which 


Peedi о пка Шене ggins (1970) is given by 


sina 
b 


f D 





av is the longshore component of The mean current velocity. 
БӘМЕШ о о а) inornron (1969) anda Longuet=Higgins (1970) 
Бит оштес The generation of longshore currents, due to an 
oblique wave approach, to The longshore component of the 
momentum flux (radiation stress) of the water waves. The 
emeceulation Of the wave-induced !ongshore current velocities 
and changes in mean water level requires the specification 
пе radiation stresses as a function от the location and 


wave properties in the nearshore region. 


Smee WAVE SET-UP INSIDE THE SURF ZONE 

As waves approach the coast and shoal there is a change 
in the momentum flux of the waves which is balanced by a 
change in thé mean water level. Outside the surf zone there 
is a set-down while inside the surf zone, after breaking, 
There is a set-up or superelevation of the water level. The 
wave set-up is important because both the local wave height 
and speed are functions of the total local water depth which 
is unknown. The change in mean water level required to 
balance the excess momentum flux of the waves must be 
determined first. A convenient form of The x-component 


(shoreward) of the momentum flux equation integrated over 


ШӘ 





ЕЕЕ ШШ еге (порите, derived by Longuet-Higgins and 


Stewart (1962), for describing the wave set-up is given by 


95 5 


> + pg(h+n) 


LS 
— 0 (1) 


Ox 





which says that the change of excess momentum flux due to 
wave action ("radiation stress") is balanced by a change in 
The mean water level. It is assumed in the derivation that 
the net local mass flux perpendicular to shore is zero so 
BU "реге (5 no contribution from The mean motion to the 
momentum flux perpendicular to the shore, and that the mean 
stresses are negligible. 

[Inside The surf zone it is assumed that the radiation 
STress Tensor can be expressed in terms of the energy and 
wave speed in the same form as in shallow water. This implies 
mien tne breaking waves are of the spilling type and that 
even under breaking waves the water particle motion retains 
much of ¡ts organized character as described by linear wave 
r Using the shallow water approximation.that the group 
Mella [Ty equals the phase velocity and that the angle of 


incidence equals zero, the radiation stress term reduces to 


"E 2 - 2 
IE = 99 Oey hah) . 22 


 Бетітитіпс (2) іпіо (|) апа іптТтеагатіпа gives the mean 


Water elevation of the form 


a KCh, 7h) + ПЕ es) 


20 





where 


Ка жено , (4) 


and the mean water elevation at breaking 15 given by 


Е үн, 
nc ESI uS (5) 
EE sa single breaker height which corresponds To the 


b 


significant wave height observed at the breaking position. 
For this study, it was assumed that the waves followed 


the breaking index, 


derived from a modified solitary wave theory by Munk (1949), 
Е оа local depth of water, О, is obtained by combining 


the local mean water super elevation and the local depth, 


DA (6) 


A sinusoidal description of the waves was used to solve 
mor the wave set-up in order to get a closed form analytical 
fer iOn апа то circumvent the difficult numerical solution 
of equation (1) required by the random-sea model. Collins 
(1972) compared the sinusoidal solution to the random sea 
Solution as shown in Figure 3. The effects of the random-sea 
model is to smooth the waves set-up curve. The magnitude of 


The difference between the two solutions is very small but 


D 





МАРИ О die LoS о лем JO OSI 1ECWOD” “ос. eu mh J 


MIA ои 50. Бы! ле вәлем оро лә 107 $] |hsoJ уез, ри Seth) | 
MI ре со. г Іс соо је“ алое цо Яс ап-475 одем JO UO рел кен р ү! 


13 “ЗЧОН54340 JONVISIO 


008 009 ООР 002 О 


`A E 


сел / 





O + N 
43 *dN-13S 3AVA 


g 
> O! 
¿| 


914 DWOJYI0UON 


обе = D 


| _ „0 

28 065 <-Н> 08/1. 
44 009 = 91 

ti 

E АШ 





mine Им тегесе can be large as the depth of water 
approaches zero. The total depth of water is important in 
prescribing the breaker point or the limits of integration 
on the Rayleigh distribution for the random-sea model. 
mene calculation of tongshore currents using The random 
sea model it is the area under The Rayleigh distribution 
ШШӘШ |5 used so that small errors on the limits generally 
cause only even smaller errors in the area. Hence, if is 
felt using sinusoidal wave descriptions to calculate wave 


Æi up is a reasonable approximation. 


MNAE FIELD INSIDE THE SURF ZONE 

A description of the wave field ¡s required їп the long- 
Shore current calculation because knowledge of it is needed 
КЕ ресі тупа the horizontal water particle velocities and 
for determining the longshore component of the radiation 
Stresses in an irregular wave field. Inside the surf zone 
the waves are unstable and the fluid motion loses some of 
Eus ordered character; but Thornton (1976) points out that 
most of the water particle motion in the body of the fluid 
is coherent with the surface and can be considered wave-in- 
а and пот furbulent, particularly for spilling type 
Штеакег5. 

Ея “а r nc aredd Rayleigh distribution as shown 
E i gure 4 is used TO give a statistical description of the 
wave field as described by Collins (1972) and Battjes (1974). 


The basic assumption (5 that at each depth a limiting breaker 


2:3 
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ВЕ ренает med which cannot be exceeded by the indi- 
vidual waves of the random field, and that those wave 
heights which in the absence of breaking would exceed the 
breaker height are reduced by breaking to the value of the 
local breaker height. That is, the energy corresponding to 
the height in excess of the local breaker height is assumed 
to be dissipated. The limiting breaker height decreases as 
depTh decreases. 

ОЕ ОО ИПА “ше Баутетан distribution, a fictitious, 
or reference, local energy per unit area, denoted as rr is 
defined. The reference energy density refers +o that energy 
density that would exist if breaking had not occurred nad 
ЕЦ јо for shoaling and retraction transformation. Battjes 


(1974) also defines reference wave heights p and their mean 





square value H МЖНШНСШ 1/5 related TO к, according TO 


The reference wave heights are assumed to be Rayleigh distri- 


buted. 
The Rayleigh wave height distribution is 
eloped at H = Ча in accordance with the assumption That The 


height of a breaking wave equals the local breaker height, 
Не, іш огсегс (9 сорта т ап approximation to the actual wave 
meant distribution. Then, the mean energy per unit area at 


с хеса point, taking account of breaking, is calculated from 


a 2 Е 
Е = = рӯн. (8) 


2 





ле маг лапсе is calculated from The pdf ot H, 


со 


не f ТЕ p(H) dH , (9) 
O 
where 
H 
S  — онь. 
p(H) = | dL 1-exp(-H*/H_*)] 
O PP =. 
zur — exp(-H “7/H А ЕШ 
5 © Г 
LN 
12 -| Н2аГ1-ехр(-н2/н 92 
О —_—_ _ 
Бо ЖОО oe (11) 
© 5 r 
A 7 ној а ит. 
IA IS Cie 


Пени ЫЕ Raylcigqh distribution implies That all waves 
Рот Ei то infinity that were previously larger than а now 
are reduced to the same height, п Therefore, the total 
probability (percent) of waves having the height Е is given 
оу 


со 


/ РЕН) dH 


H 
S 





Mie contribution то the variance is given by 


со 2 


Н 
ај зн АН = Н ^ ехр(- >, 
5 © H 


H r 

S 
КЕЕ 5 The Term on The right of equation (Il). ki is the 
local breaker height, which inside the surf zone is assumed 


to be given by 


П-р МЕ 


Mme local mean wave height R can be expressed in a similar 
manner in terms of ШЕ and pic by means of the clipped Rayleigh 
Шістгіритіоп, 

е E ТТІ 2 D 

Н = | Hd[ | -—exp ( -H Jin ) | + A exp(-H, шч ) 

О (14) 

where again the term on the right represents the percent of 
waves greater than F: I ено па та Јани грба тог which now 


have the height 5/27 Integrating (14) gives 


- yn 
Я = — 


Н EL A ) MESA 
r se 
rms rms 


in which 


NO 


Ке 


ae 
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and The error function being defined as 


p 
erfíp) = ыы Sume es : 
yn 
О 


ШОШ ег оГ function was calculated using the rational approxi- 


ПОО ОП ot Abramowitz and Stegun (1965), 





erf(p) = I-[(a,t + azt + Po ие, 
where 
+70), а, = n022; 
2 = .47047, a) = 025957985 
Ке нис, n 


Цане сте j terror cis ina This approximation is 


qeu ox пате 


In the observations used for comparison with The Theory, 
the breaker height is measured visually. It is assumed 
That an observer visually measures the significant breaking 
wave height defined as the average of the highest one-Third 
И О Т оп of the wave heights. Тһе difficulty in applying 
this definition to the present problem is that the defini- 
Tion applies To a point measurement or a statistically homo- 
geneous (spatial) wave field and in this problem the waves 


are defined as varying spatially as they shoal shoreward. 
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In order to define the significant wave height for a 
Spatially varying (nonhomogeneous) wave field, if is 
assumed that the observer measures waves when spatially 
one-third of the waves have broken; hence, the reference 
Mave height can be specified from the clipped Rayleigh dis- 


ШЕ 9017 | оп 
ехр (-Н үне me 5707355355 


zZ те 


a H . (16) 


Bee LONGSHORE CURRENT VELOCITY 
Леа ТИ аып от the longshore current starts with the 


y-momentum flux equation 


gor 
9х Б ан I 





ee 


Where the lateral shear stress, T is neglected (17) reduces 


es 
ШО 


95 


ху E 
3x t Шо 0, (18) 





ШЕП саус Тпат the change of y-component (longshore? 
momentum flux due То waves in the x-direction is balanced by 


the bottom stress, T сш = ест оп, Assuming That 


b7 


the amplitude of the wave motion |Uw| is much greater than 


Mmmemmean current velocity, V, then (Thornton, 1969) 


E Ud (19) 
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The excess of momentum flux of the waves, or "radiation 


stress" component is given Бу 


E 
S =E PS [nG соба 28 
ху 2 


, 


which inside The surf zone reduces to 


S = E singa cosg (20) 

ху 
under the assumption that Со = C for shallow water. Combin- 
ing equation (8) and (12) gives the mean energy per unit 


area as 


p 5 ра H БОЛСОК А 2/91 (21) 


and 


m —— | 
r= 


ра H m ss E осо: „е 0223 


| 
> xy 8 


Тһе variables "a and à can be expressed as differenriable 
ml ions of x (distance from shore). Recalling that it was 


assumed inside the surf zone 


сүр, 


muere the total depth is the sum of the still water depth 


plus the set-up 


E 





oie ata (ОП (5), (5), (6), and (13) ДР can be expressed 


as 
y? 
ba > H, (K - Tg) E ce) 023) 
Application of Snell's law allows the local breaker 


mole, a, tO be expressed in terms of the known breaking 


alle а апа the breaking celerity C 


b D: 


NT. (24) 


Using the shallow water approximation for wave speed 


n= e 
and 
b g Y 4 
Then 
а (25) 
EN b 


memce, trom equations (23) and (25) it is seen That both 
ive aac Quer mow Expressed as functions of h, which in turn 
EN function of x. 


The bed shear stress coefficient is determined by com- 


Ила едџат!оп (18) апа (19), 


18, | 
5. - > 


Fr a (26) 
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lu теал погодат water particle velocity amplitude 





is 
expressed using linear theory (Battjes, 1974), 
Џи] = ==, (27) 
where H is given by equation (15). 
The change in the radiation stress is given by 
ER e 
95 
EN s ez. bs 
Ox Zl b H H 
b b b 
+ A SO exp ( -H 21H 2 028) 
4 5 5 r à 
where 


3 x 
ЕЕЕ i [| ехр(-Н. 


2 2 
/H. 


) ] 


and 


vos ads 
sina cosa = ¿ME e sina, (h - ER Y 
Hp b H, 
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O SATA 


ШЕ LEO DATA 

Darron ем ико па Environment Observation (LEO) pro- 
Gram established by the Coastal Engineering Research Center 
was used in this study. In the LEO program, nearly simul- 
Taneous observations of breaker conditions (height, period, 
mere Of approach and Type), local winds, ltongshore currents, 
foreshore slope, width of the surf zone and rip currents were 
made daily during fhe period under consideration. The long- 
Shore current was determined by observing the direction and 
measuring the distance parallel to shore thar a dye packet 
injected into the surf zone traveled in one minute. Appendix 
A provides the set of instructions followed during the obser- 
i ions. 

чата used тог this study cover a period from May 1972 
e ember 1975 and refer to stations: 5705, 5706, 5707, 
ЕЕЕ 2714 and 5715, located within the confines of Point 
Mugu Naval Air Station, 60 miles northwest of Los Angeles, 
Niro nia (location 6 їп Figure 5). 

From these stations, 4,6532 observations were considered 
of which only 647 data points were used in the analysis. The 
MONO wing criteria were discussed to eliminate observations 
which were not consistent with the application of the theory 


Or were simply erroneous: 
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Rip Currents 

lt was beyond the scope of this study to account for 
o ie tec ог The longshore current system by rip 
 Іггейпісв. Hence, all observations noting the presence of 
EM ШШ ел еге systematically deleted. 

2. Angle of Wave Approach 

The theory used in The derivation of the formulas 
Fuss im this study assumes that the angle between The 
direction of wave approach and the depth contour must be 
different from 90 degrees in order for a longshore current 
то be generated. Hence, all observations іп which the wave 
direction was reported as being perpendicular to the shore- 
line were neglected. 


Di, Wind 





Shepard and Inman (1950) suggested the importance 
Пе wind Im generating longshore current; They also indi- 
Gate That it is difficult to separate The wind generated 
current contribution from the current generated by The waves. 
Thus, observations where the wind speed was reported as 
being greater than ten miles per hours were not considered. 
4. Foreshore Slope 
Observations where the foreshore slope was reported 
as being greater than ten degrees were neglected since such 
large values are not consistent with what is usually observed 
ЕШ (Пе beaches under consideration. 
Ec Wave Period 
аи O кебер туе study restricted fo sea and 


swell of relatively short period, all observations where the 
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Wave period was reported as greater than 20 seconds were 
neglected. 
G Done Бага 

All observations in which the reported data were 
considered to be incorrect due to either mistakes of Тһе 
observer or the typist, such as l|ongshore currents in ex- 
Mes Of six feet per second, direction of approach greater 
than 180 degrees, distance of dye injection greater Than 


ЕШ) Теет, etc., were systematically rejected. 


КОИ OUIRCES OF ERROR 

(Рошела ine winnweres: in longshore currents, if is 
somewhat surprising that there are relatively few sets of 
adequate field measurements of longshore currents and the 
simultaneous wave parameters in The surf zone. After a 
Search of тһе literature it was concluded that little has 
been achieved for devising electronic equipment designed 
for gathering longshore current and associated wave informa- 
п оп а routine basis. Hence, as in The case of the LEO 
ES most of the observations must rely on the good 
judgment and personal abilities of the observers. This 
luEOduces a subjectivitfy factor which ultimately affects 
me tinal results. 

|. Breaker Angle 

Galvin and Nelson (1967) suggested that the variable 

most difficult to measure with necessary accuracy is the 
emgle of wave approach or wave direction. Galvin and Savage 


(1966) suggested that when using a visual compass referenced 
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to a baseline to measure the breaker angle, the errors may 
ОНИ Бем nwo degrees, leading to a relative error which 
is very large for smal! breaker angles but which decreases 
as the breaker angle increases. 

In the LEO observations a protractor was used for 
determining the breaker angle as shown in Appendix A. This 
system is completely visual using the unaided eye to estab- 
lish the perpendicular to the shore and introduces a greater 
mumen factor. Hence, if is a good assumption to attach an 
accuracy less than that suggested by Galvin and Savage to 
Such measurements. 

en Beach Slope and Surf Zone WidTh 

In describing longshore currents which flow within 
The surf zone, accurate knowledge is required of the beach 
Hume including both the beach slope and the width of the 
surf zone. The approach used for the LEO data was to assume 
a plane beach inside the breaker line. The nearshore sub- 
aqueous slope was computed using the observed surf zone 
width and the observed breaker height. Thus an uncertainty 
Numer for The beach was introduced. Referring to LEO obser- 
n instructions in Appendix A, the observation of the 
ee ZOne width "... is based upon the judgment of the obser- 
Iu Man-made or natural features in the surf (e.g., a pier) 
да in This observation." Again, a subjectivity factor 
involved. 

The observed foreshore slope could not be used because 


it proved TO be unrepresentative of the beach profile within 
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the surf zone. The calculated water depths inside the break- 
пе Using This Slope were systematically greater than 
the calculated breaker depths. 
3. Wave Period 

Galvin and Nelson (1967) point out that under favor- 
able conditions the wave period can be measured with reason- 
able consistency in the field by visual observation. Although 
mms Parameter was not used directly in the computation if 
is interesting to notice that their suggestion agrees quite 
well with the LEO observations since the range of periods 
found fall into the expected values for the shoreline under 
sonsideration. 

4. Breaker Height 

In the LEO program, the breaker height observation 
is based solely on the judgment of the observer. Known 
dimensions of natural or man-made features on the shoreline 
or in the surf zone are used as references for estimating 
the wave breaking height. Galvin and Savage (1966) sug- 
gested a relative error in breaker height measurement of 
NE о – гсепе They arrived at this figure by comparison of 
breaker height measurements made with pressure gages, 
oscillographs and visual observations, although the measure- 
ments were not made simultaneously. Hence, in the light 
of their finding it can be concluded that at least the same 
Eu should be expected in the LEO data in which the obser- 


Ex ons are solely visual. 
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Би опао пеге Current 
In the LEO program, the current speed was determined 
HEN Козы е ае а tracking agent. This also adds an uncer- 
uv factor due то fhe diffusivity characteristic of the 


dye. 


07 СВЕ ВАТА 

A set of 62 field observations made along the Southern 
Ser trornia coast (Fig. 5), taken from Ingle (1966) and per- 
sonal communication, were selected using the same criteria 
used for selecting the LEO data. Despite the size of Ingle 
Sample, about теп percent the size of the LEO data, ifs 
analysis is important since the Ingle observations are more 
accurate Than the LEO observations. Thus, the Ingle results 
serve as a reference comparison to the results obtained using 
me, LEO data. 


The parameters H period, a, and V were taken directly 


b b 


from a summary appendix in Ingle (1966). The beach profiles 
апа Тһе distance shoreward from breaker zone in which The 
longshore current velocities were recorded also were avail- 
able. However, the positions of the breaker were not avail- 
able. Thus, the parameter np was obtained from the relation- 
Ship Jm = 1.28 Hp; and the local depth below still water 
level, h, and the beach slope were scaled out from the beach 
profiles presented in the publication. 

It should be pointed out that in the Ingle data the break- 
er heights were measured by sighting on either a graduated 


Ее held at an approximate still-water line and the horizon, 
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a graduated pole held in the zone of breaking waves, or a 
piece of cardboard with a slit and graph paper along one 
edge. Breakers less than 2 feet in height were estimated 
while standing in the breaker zone. о E ONDES 
breaker height observations were based solely on the judg- 
ment of the observer on the shore. For measuring the 
breaker angle, Ingle observers used a Brunton compass while 
Standing in the surf zone, supplemented by sights taken 
from positions elevated above the beach; LEO observers used 
o racior, aS Shown in Appendix A, with the observer on 
шее beach. For measuring other parameters, both LEO and 
Ingle observers used essentially the same techniques. аг 
КОИ трогтапт to mention, that in the case of the Ingle 
Observations, the beach slope and the position of the obser- 
vations in the surf zone were better than those of the LEO 
observations since in the former an ordinate and abscissa 
arrangement of.wooden stakes allowed workers to position 
Themselves in the surf zone; also most people involved in 


Ingle observations were well trained personnel. 
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SE SUL TS 


БЕС 222 (907 (ој, 512), (15), (16), (21), (22), 
БЕ 25) (260), (27) and (28) presented in the theory sec- 


tion were used with the parameters H x and У from 


PCS 
LEO and Ingle field observations to solve for the bed shear 
sess Coefficient. The coefficient calculated is based on 
Eu specific locations within the surf zone and not for 
men conditions. It should be mentioned that since the co- 
efficient is not determined by direct measurement, itt there- 
Meme пот Only reflects bed shear stress, but also any errors 
and uncertainties in measurement. 

Ба е the assumption that The observer 
measured waves correspond То а clipped Rayleigh distribution 
when spatially one-third of the waves have broken, used for 
computing the reference wave height Низ was tested using 
other assumprions. тера пер ас logically be argued 


that the significant wave height might correspond to the 


point where half the waves have broken, in which case 


一 一 一 一 一 


2 Zu = 
exp(-H, Zur 250/155 


This assumption was used for computing new values for the 
coefficient. The relative difference between the coeffi- 
cients Thus calculated and the original ones was determined; 


ШЕР Уагізсіііту мас Тоипа To be of the order of ten percent, 
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which is relatively low. Hence, it can be said That the 
bed shear coefficient calculations are not very sensitive 


to the assumed definition of the significant wave height. 


КОО LEO DATA 

Figure 6 depicts the frequency distribution of values 
obtained for the coefficient and selected statistics of the 
ет“ граот топ. The variability in coefficient values, as 
represented by the standard deviation of the distribution, 
Na measure of the consistency of the calculation of the 
coefficient from the field observations. The mean of the 
distribution is 0.008, while the standard deviation is 0.010. 
This suggests there is a large spreading in the results. 
Be cver. iT should be noticed That more Than 90 percent of 
The calculated values fall between 0.001 and 0.020 and that 
Nena istTribution has less spread than a Gaussian distribu- 


tion for fhe same standard deviation. 


Er INGLE DATA 

Figure 7 depicts the frequency distribution of values 
ишы леа for the coefficient and selected statistics of the 
Meer ibuUTIOnN. The mean of the distribution is 0.014 and 
the standard deviation 0.01], which suggests again a large 
spreading in the results. However, 95 percent of values 
EL Nbeiween 0.001 апа 0.030 and again the distribution has 
less spread than a Gaussian distribution for the same 


Standard deviation. 


42 





200 


180 


160 


140 


к 
O 


со 
O 


> 
© 
A A НЕЙ 


FREQUENCY 
© 
© 


60 


~ 
© 


Br. 


ПО 015 0921 026 031 036 040 045 050 055 060 


COEFFICIENT 


Mean 0.0908 Skewness 5.694 
Yarlance 000009 Kurtosis НО 
га. Vevrat ion 0200972 Minimum 2008 
бое, га Гоп ЕЭ Maximum Ore Oe 
Figure 6. Frequency Distribution of Coefficient Values 


ОЕЕО рата). 


43 





12 
11 
10 
: | 
> 
8 
5 | 
= 
ul 7 
ES 
© 
u 6 
rr 
ш 5 
4 
3 
2 
| + 
‚001 008 016 1O23 O31 038 045 
COEFFICIENT 
Mean OO Skewness 1984 
Variance (ees Kurtosis 2.649 
ова и кај Гор OSO 2 Minimum OS 001 
(ере о о“ a Lon OS OS Maximum OROSO 


MUS 7. Frequency Distribution of Coefficient Values 
(Inale Data). 





COMPARISON OF RESULTS 

Despite the difference in sample size between the LEO 
and Ingle data used for the calculations, some comparisons 
can be made. A simple way of comparing both sample results 
is by looking at their mean and standard deviation. The 
mean and standard deviation corresponding to the values of 
The coefficient for the Ingle data are both larger than the 
values obtained for the LEO data. There is a relative dif- 
ference of /5 percent between the mean of the coefficient 
values of the Two samples; but, the relative difference 
between the two standard deviations is only ten percent. 
This says that the distribution values for the coefficient 
In both samples is nearly the same, although for The Ingle 
data the values for the coefficient were somewhat larger 
and with more spread than for LEO data which might be ex- 
pected for the smaller sample size. 

БО ае ыле when comparing bor sets of data the 
Ingle observations were more accurately taken and more 
reliable Than the LEO observation. Hence, the results ob- 
Tained with the Ingle dafta would be expected То be better 
urhe results obtained with the LEO data. То test if 
there is any statistical difference between the two sets 
of data, a hypothesis test about the two means obtained was 
ED тле central limit theorem states that, if x is the 
team OF a random sample of size n Taken from a population 


пам па rhe mean U and the finite variance 9°, then 
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Ден а |ш= ог а гепдот variable whose distribution function 
Васо лез Gnas Of ће standard normal distribution as n > e, 
The variances of the population are unknown, but since both 
samples are fairly large, it is justifiable To approximate 
ви етой variances with rhe samples variance. Thus, a 


Æ i Statistic can be stated as 


The hypothesis to be tested is the null hypothesis, Hu, = ОЕ 
against the alternative hypothesis р-у > 0, where ы герге- 
sents the mean of the population. The evaluation, for the 


data available, of the z statistic was found to be equal to 


EN 00S + = 42-525 
(9.0001 - 0.00009, 1/2 Ë Y 
62 | Ba 


Ma level of се и ООС ОТ ОТ ОО! The z statistic for The 
Папа distribution is 5.49. Since the value obtained for 
ет statistic is larger than the critical value of 3.49, 
The null hypothesis is rejected with great confidence; and if 
can be concluded that the difference between both means is 
st ically staniticant and cannot be attributed to chance. 
Therefore the results obtained with Ingle data are better 


шпап rhe results obtained with the LEO data. 
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ICO RRE TION МАТЕ INDEPENDENT VARIABLES 

ATTempis were made To Correlate the calculated coeffi- 
cients with the independent variables, breaker type and wave 
period, which were recorded in the field but which were not 
used directly in the computations. Analysis showed nothing 
Sonciusive regarding The correlation of The coefficient to 


Me breaker Type since The distribution of breaker types 


among the data was very uneven; the spill/plunge type repre- 
sented 7/2 percent of the data and the spilling type 20 per- 
cent. Table Il shows some selected statistics of the ais- 


tribution of coefficient values for various breaker types. 


IEEE clecsegoSTaristics-ftor Distribution of 
Coefficient According to Breaker Type 


(ево Waid ) 

No. O? as ina surging “sopill!/Plunge 
Observation (50 27 25 454 
Mean O50090 000,5 0.0098 000738 
Variance ОШООО! 0.00004 07000155 0.9009 
ENS Dev. 0. 0107 0.00609 0.01204 0.00946 
Met. Var. ево OO 9 On 52 Река Зајко 
Range 0.0790 000250 0205500 C207 200 
Minimum оого 0.09 70 0.001200 0.00100 
Maximum 020200 0.0240 0.05400 0075800 
Skewness x20 БЕС! 209024 506 О 
Kurtosis 176.0 10555 7.38989 во: 
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A simple linear regression between the calculated coef- 


ficients and the observed period gave the selected statistics 


Бг Table |||. 
TABLE ШЕКЕ СОГ Г ella on of Coefficient with 
Wave Period Statistics 
EEO бата ОЕ Юта 

Correlation CR) =E 03570 - .08479 
Sid. error of estimate . 00 о, aol ee 
R squared -00290 -00719 
ошаптт гсапсе . 05574 ООС] 
Intercept 20059 0.892 
Slope = коор, - „00042 


The negative sign of the correlation coefficient indi- 
cates that there is an inverse relationship; that is, the 
value of the coefficient tends to become smaller as the 
period increases. However, this relationship is very weak 
dindi cated by The absolute value of R which in both cases 
Exmuch smaller than one. This result is not surprising 
since waves in shallow water become non-dispersive or in- 


Шо апт of period. 
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YRC ONC LUS TONS 


An analytical solution for the bed shear stress coeffi- 
cient was derived using the concept of radiation stress. 

The best theory for calculating the variation of wave 
energy and longshore current, and the resulting bed shear 
Bes coefficient, was to use The truncated Rayleigh p.d.f. 
for the statistical description of the wave field inside The 
surf zone. A sinusoidal approximation of the waves was used 
Киса соате the wave set-up. Calculations of the coeffi- 
cient were made by using suitable sets of data obtained 
ла the LEO observation program and Ingle (1966) observa- 
шета along The Southern California coast. 

ет ари ту im he results obrfained for the coefficient 
values were expected due to subjectivity and uncertainties 
in the techniques used in The data collection. This ts The 
MMS test of the bed shear stress coefficient using fairly 
large setsof field measurements within the surf zone. Even 
with the uncertainties involved, the analysis resulted in a 
fairly good agreement between the mean of the calculated 
coefficients in this work and the values obtained by various 
investigators for the bed shear stress coefficient for dif- 
ferent test conditions and outside the surf Zone. 

|t was shown that the dependence of The coefficient on 
the wave period is negligible, in agreement with the assump- 


tion that waves inside the surf zone are non-dispersive ог 
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ШОШ ОООО Л аш Since one of The biggest differences 
between Pacific and Atlantic coast waves is the period, if 
feb eee Concluded That fhe calculated coefficient is not 
Ocean dependent. 

Since it was initially concluded that Ingle's data was 
Of higher quality than the LEO data, it is assumed The 
coefficient values using Ingle's data is therefore more 
reliable. In any event, the mean value of the two data 
sets are the same to two significant decimal places. There- 
fore, it is concluded that a reasonable value for the bed 


МЕНИ (гео сое: стеле within the surf zone is 0.01. 
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APPENDIX A 
LITTORAL ENVIRONMENT OBSERVATIONS 


CERC Farm No. 113-72-8 Mar 72 has been designed for keypunching onto computer cards directly (small numbers above each box represent card colum numbers). 
It is recommended a pencil be used. All data should be recorded carefully and legibly. Errors should be corrected by first erasing erroneous data as 
write overs usually produce illegible data. Make remarks as necessary on the form but record only data ir, the boxes provided. All observations must be 
made at the same point on the beach every time (in front of the reference pole). 


STATION IDENTIFICATION: 
Each site in the “Littoral Environment Observation" study has been assigned a numerical code consisting of 5 digits. The first wo digits define the state 


or territory in which the site is located and the remaining 3 digits define the particular beach or park within the state or territory. A space has also 
been provided to write in the name of the particular beach or park at which the observation is taking place. 


DATE: 
Indicate in the spaces provided the year, month and day on which an abservation is made. 


TIME: 
Indicate the time at which the observations are being made. The 24-hour system of recording time has been selected in order to eliminate any confusion be- 
tween AM and PM. The hour "OO" refers to midnight, "O7" to 7:00 AM, "13" to 1:00 PM, etc. 


SURF OBSERVATIONS : 

a. Wave Period - Record the time in seconds for eleven (11) wave "crests" to pass same stationary point. Eleven "crests" will include ten camplete 
waves (crests and troughs). The first (1) "Crest" selected for observation is recorded as time zero and the eleventh (11) "crest" will be the stop or cut 
time. Record this time in seconds in the spaces provided. 


b. Breaking Wave Height - This observation is based solely on the judgment of the observer. Natural or manmade features om the shoreline or in the 
surf zone whose dimensions are known may aid in judging the height of a wave. Otherwise the cbserver's best estimate will be sufficient. Record the breaker 
height to the nearest tenth of a foot. 


c. Breaker Angle - To determine the direction fram which the waves are approaching the beach use the protractor on this reverse side of the data form. 
The 0-180? line should be oriented along the shoreline; use the protractor to site the direction from which the waves are approaching when they are first 


breaking. 
d. Type of Breaking Wave: 


Spilling - Spilling occurs when the wave crests becomes unstable at the top and the crest flows down the front face of the wave producing an irregu- 
lar, foamy water surface. (see figure 1) 
Plunging - Plunging occurs when the wave crest curls over the front face of the wave and falls into the base of the wave producing a high splash and 


much foam (figure 2) 
Surging - Surging occurs when the wave crests remains unbroken while the base of the frant face of the wave adcances up the beach -ftsee figure” 3r . 
Spill/Plunae - A cambination of both spilling and plunging occurring simultaneously. 


WIND OBSERVATIONS: 


а. Wind Speed - A wind meter is provided to each observer and it is recamnmended that the instructions orovided with the meter be followed to obtain 
wind speed measurements. 


b. Wind Direction - After the approximate orientation of the beach with respect to north has been defined the observer can determine the direction 
“fram wnich" the wind is coming. 


FORESHORE SLOPE: 

For measurement of the foreshore one must use either the clipboard/inclinameter or the Abney hand level. Cbservatians should be made as close to mid-swash 
as possible. Using the climboard/inclincmeter place ıt on the appropriate edge and record the angle where the ball cames to rest. Using the Abney hand 
level place it on a straight edge and level the bubble; record the indicated angle. 


WIDIH CF SURF ZONE: 
This observation is based solely on the judgment of the observer. Estimate іп feet the distance from the shoreline to the line of the mst seaward breakers 
(not to be confused with white caps). 


LONGSHORE CURRENT: 


a. Dye Distance - Dye packet should be injected just shoreward of the breakers, if possible. Driftwood or any other floating object should be used if 
dye is not available. Estimate the distance from the shoreline to point of injection and record this distance in feet. 


b. Current Speed - Mark the beach in line with the injected dye and make a second marx to indicate the dye movement after one minute has lapsed. Pace 
the distance between these marks and record this distance in feet. 


с. Current Direction - When looxing seaward, if the dye has moved to the left record -1, to the right record +1, and no langshore moverent record 0. 


RIP CURRENTS: 

Rip currents are defined as seaward moving channels of water which return the water that has been piled up along the shore by incoming waves. Rip currents 
are fed by feeder currents, water moving along the shore (see figure 4). Two currents join and extend out in what is known as the "neck", where the water 
rushes through the breaker zone in a narrow lane. Beyond the breaker zone the current spreads out and dissipates in what is called the "head". If such 
rip currents are present estimate their spacing in feet. If no rips are present record 0. 


BEACH CUSPS: 
Cusps are semicircular or crescent shaped cutouts in the beach face (see figure 5). If such shapes are observed record the distance between the "hams" of 
the cusps which indicate the spacing. where the spacing is irregular estimate the average spacing. If no cusps are present record 0. 


FIGURE 2. 
PLUNGING WAVE 
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